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Abstract 

The equivalence of the Rivier-Margenau-Hill and Born- Jordan- Shankar a phase 
space formalisms to the conventional operator approach of quantum mechan- 
ics is demonstrated. It is shown that in spite of the presence of singular 
kernels the mappings relating phase space functions and operators back and 
forth are possible. 
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1 Introduction 



A wide class of quantum quasi-probability distributions F of position q and 
momentum p was studied by Cohen the Wigner function and the 
Margeneau-Hill functions || being particular cases. The choice among the 
F functions representing the quantum state of a system is similar to the 
choice of a convenient set of coordinates. However in order to use safely one 
of these phase space functions it is necessary to show that the usual operator 
formulation of quantum mechanics is equivalent to the corresponding phase 
space formalism. In the later, in addition to the state, it is necessary to 
specify the functional form of the dynamical variables. In Cohen's approach 
each of the distributions F is obtained with a particular kernel function / 
from the quantum density operator g, g — > F. (Operators are represented 
with an accent and the notation corresponds, as in P], to a particle 
moving in one spatial dimension.) Cohen gave as well the quantization rule 
g — > G that associates with a classical function g(q,p) a quantum operator 
G(q,p) in such a way that the expectation value of the operator can be 
written equivalently as a trace or as a phase space integral (All integrals in 
this work go from — oo to oo.) 

(G(q,p))=ti(gG)= J! F(q,p)g(q,p)dqdp. (1) 

The most frequently used quantization or ordering rules and many others 
fit into this general scheme (in particular the rules by Weyl [§]||, Rivier jj, 
Born- Jordan J/J, and the set of rules known as normal, antinormal, standard 
and antistandard P-|TD|), together with their associated phase space quasi- 
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probability distributions, see Table 1. 

In order to build a quantum formalism in phase space equivalent to the 
conventional operator approach, rather than considering g and g as primary, 
a different point of view is required where the primary objects are the state 
and observable in operator form, £ and G. Their phase space representations, 
F and g, are obtained from them using g —>■ F and the inverse transformation 
of the quantization rule, i.e., G —* g. Also in this case it is imposed that the 
expectation value of G(q,p) is given by (P, but now the function g(q,p) is 
not necessarily equal to the classical magnitude; it is simply one the images 
or representations of G in phase space. 

It is also possible to consider F as a primary object. Then, the trans- 
formation F — > g is required to obtain the corresponding density operator. 
Even though it is not common practice to consider F as primary in quantum 
mechanics, there are physical systems of practical importance (semiconduc- 
tor heterostructures) which are modelled in this fashion flIT| , |T2|| . It is also of 
interest to note that classical statistical mechanics is a theory formulated in 
terms of an F distribution (the classical distribution function) and proper- 
ties g (classical magnitudes), so that the transformations F —>■ g and g — ► G 
provide a set of equivalent operator formulations of classical statistical me- 
chanics. Their potential applications are yet to be fully explored. One of 
them, using f — 1 and associated with the Weyl-Wigner formalism, was 
examined in [p~3|j . 

The (formal) full set of transformations that completes Cohen's two orig- 
inal mappings, has been described by several authors |14]-[nj. Whereas the 
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transformations g — > G and g — > F involve the kernel function / in a multi 
pie integral, the inverse transformations involve / . Agarwal and Wolf WE 



restricted their detailed study to mappings where the kernel function had no 
zeroes to avoid the singularities of f~ x . In fact it has been generally believed 
that the inverse mappings cannot be performed in these cases ||17|| . As a 
consequence the investigation or applications of some of the / functions and 
their associated quasi-probability distributions and ordering rules have been 
scarce. Our main objective here is to demonstrate that these singular kernels 
do not necessarily preclude the existence of the inverse mappings and that 
indeed equivalent phase space formalisms based on them can be constructed. 
In other words, we shall show that there is an "inverse operator basis" , see 
(15) below, in the form of operator valued distributions, associated with 
these kernels. The need to consider in general the mappings between opera- 
tors and phase space from a generalized function point of view was already 



emphasized by Agarwal and Wolf [15 



47T 2 J J J \ 2 



u - ll) e -mq-u)+r P ]^ e ^ r ) d9 dr du (2) 



The quasi-probability distributions in phase space are obtained from g as 

Q 

Cohen noted that by imposing the condition 

/(O,t) = /(0,O) = 1, (3) 

the resulting F function provides the correct "marginal distributions" for 
q and p. A group of / functions and their corresponding quasi-probability 
distributions are listed in Table 1 |Tj|. The property (||) is desirable but 
this condition is not fulfilled by several useful quasi-probability distributions 



[such as the P-distribution by Sudarshan || and Glauber ]B| or the Q (or 
Husimi |P0[| ) distribution.] 

The operator G(q,p) is given from the phase space function by 

G(q,p) = ^JJJJ 9(q,p)f(e,T)e- l ^ + ^-^dqdpd6dT. (4) 



It is an exercise of Fourier transforms to obtain g in terms of G from 

G 



g{q,p) = 7T f f f ( u + 



2ttJJJ \ 2 
Similarly 



rh \ e~ i ft( q ~ u ^ +Tp } 
U ~Y/ f(-9 -r) d9drdu ^ 



q=1LJJJJ F(q } p)f(-9, - r )- 1 e-W*-* )+T < p -d] dq dp d6 dr. (6) 

The explicit expressions for all the transformations in equations (0), (|j), (|^) 
and (H) can be summarized as 

F = A[f]g g={A[f}}- 1 F (7) 

g = hA[f- l ]G G = h- l {A[f- 1 ]}- 1 g (8) 

where / = f(—9, — r). 

Here we consider the transformation from g to F as the reference mapping 
represented by the "operator" (on the space of density operators) A[/] which 
depends functionally on /, see (@). The inverse operator {A[/]} _1 acts on F 
to provide the density operator g. Note that G — > g involves, except for a 
constant, the same operation as g — > F but a different kernel, namely f^ 1 
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These are of course formal results and for every / it is necessary to study 
if these integrals exist and to determine domains where the transformations 



can be performed. Seemingly functions having zeroes, such as f(8, r) = 
cos(6rh/2) or 2 sin(#T^/2) / '9rh, may be problematic because of the presence 
of the inverse of / in the integrands of the transformations G — > g and F — > g. 
It is also possible to relate operators and phase space using a framework 



complementary to Cohen's [ 14 -|17[ : Assume that there is an operator basis 



B(q,p; q,p) such that the operator G(q,p) can be given as 

G = If 9B(q,p)B(q,p) dqdp, (9) 

where the "coefficients", g B {q,p), 3X6 , ctS before, the transform, image or 
representation of the operator in that basis. For a basis B to be practical it 



must have an inverse B such that [^2 



tr[B-{q,p)B(q',p')} = 5(q - q')5(p - p'). (10) 
If the density operator is expanded in the inverse basis, 

g = J J F B (q,p)B-(q,p)dqdp, (11) 
with expansion coefficients FB{q,p), using (0) it follows that 

ti(Gd)= [[g B F B dqdp (12) 



and the coefficients (phase space representations of the state and the observ- 
able) are obtained from the operators by taking the traces 

F B (q,p) = tr(£B) (13) 
g B (q,p) = tr(GB-) (14) 
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In terms of / these bases are formally given by 

B-(q,p) = e-m^+rip-P)] J dedT (15) 

B(q,p) = -^JJ e- i[e(q ~^ +T{p -P )] f(9,T)d9dr. (16) 

From this perspective each of the phase space formalisms is based on a certain 
"coordinate system" or "basis" that can be more or less convenient depend- 
ing on the case. But, as before, the presence of f~ l in the expression for 
the inverse basis seems to impose a limitation when / has zeroes [ Hp . In 



this way one could be prematurely tempted to discard, for example, phase 
space formalisms based on the relatively common Rivier quantization rule 
(or symmetrization rule) and the associated Margeneau-Hill function or on 
the Born- Jordan quantization rule and the associated Shankara distribution 
function ||23|| . However the next section shows that the zeroes of / are not 
actually a problem. 

2 Effect of zeroes of r) 

We shall give an explicit expression of the transform of (pp™ 1 for arbitrary 
nonnegative integer values of n and m. According to equation (|5]) the /-image 
of (pfj" 1 is obtained by the integral 

h r f f I rh 



g(q,p) = 9(w,TiT;f) = ^ /// ( u 
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rh 

U + Y 



x e i[0(q-u)-Hp]_J_ dedTd ^ ( 17 ) 

J [yi T ) 

Introducing a closure relation in momentum this takes the form 

= ^JfJJ [u-^yp^e^-^e-^'-^jy^dedrdudp'. (18) 



Making use of the integral expression of the Dirac delta and of its m-th order 
derivative 8^ the integral in p' is solved, 



g(Q,p) 



tTi \ n 1 
'/ - — I e i0(q ~ u) e iTp — ; 5 {m) (r) d9 dr du. (19) 



f(0,r) 



The r- integral is carried out next. It is necessary to consider the m-th deriva- 
tive of the function 



2 J 



f(0,r) 



(20) 



By using Leibniz's formula for the m-th derivative of a product twice and 
putting r = one obtains 



d m X 
dr m 



E 

1=0 



m 
I 



i 

E 

3=0 



(ipY 



111 



(n + j-l)\ 



u 



n-l+j 



(-h/2) 1 



_j d m ~ l f - 1 



dr 



m—l 



T=0 

(21) 

which is used to integrate over r with the aid of the derivatives of the delta 
function. Then the w-integral is carried out, 

m I i 



g{q,p) 



:-i) n EE 



m 



I 



nl 



:(V2) 



l-j 



Jm—l f—l 
iBq a J 



dr 



in— l 



r=0 



5 {n - l+J \9) d9. 



(22) 



The derivatives of f^ 1 with respect to r at r = can be performed by using 
the formula for the derivative of arbitrary order, say m — /, of the inverse of 
a function [24]. This formula involves / _1 (t = 0), 

d™- 1 !- 1 , l 



r __ r(m - m - ir - fm _ i+i(T = o) 



D 



m—l 



(23) 
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where 



D m -i = 



a i 
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(24) 



and 



0"n = o, n {9) 



1 d n f 



n\ dr n 



T=0 



(25) 



Finally the integration over 9 is performed using the derivatives of the delta 
function, and taking the condition (j^) into account. Again, Leibniz's theorem 
is used to arrive at the lengthy but explicit expression 



m I n—l+j 

(-l) m_i {m-l)\i j ~ k - m (h/2) l - j 

1=0 j=0 k=0 

rl k D 

. pign-l+j-k u ^m-l 



m 
I 



n — l+j 
k 



TV. 



d9 k 



) (n + j-l)\ 

where only mixed partial derivatives of the / function at the origin (9 = 
0, r = 0) are to be considered. 

We conclude that provided that Cohen's condition @ is satisfied and 
/ is analytical at the origin g{q,p\q n f) m ] f) is well defined in spite of the 
possible zeroes of /. The expression corresponding to the operators in reverse 
order, g(q,p;p m q n ; f), is the same except for a the change of sign — > 
(-ft/2)K 



(26) 



Examples of images of q p for several / using 



are 



h 2 

g(q,p; fp 2 ; f = 1) = p 2 q 2 + 2ihpq - — 
q,p; q 2 p 2 ; f = cos(9rh/2) = p 2 q 2 + 2ihpq 
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q,p;?p i ;f = 2sm(9Th/2)/(9Th) 
q,p;fp*;f = e^m 



ti 2 

2 2 r\ • 4- 

p q + Ziripq 

3 

2 2 

p s 



pV + Aihpq - 2U 2 (27) 



It can be checked -via eq. (f|) with the corresponding / for each case- that the 
operator obtained from these phase space functions is indeed q 2 p 2 , and that, 
once / has been chosen, the relation between operators and phase space 
images is biunivocal for arbitrary (non-negative) values of m and n. It is 
possible in summary to map into phase space functions at least operators of 
p and q in polynomial form or given by expansions in (pp™ or p" 1 ^ 1 . A broad 
set of bounded operators admits such expansions |I(J • 

Note that ( p6|) is non- linear in / so that if f — afi + 6/2, m general 
<7/ 7^ agf x + bgf 2 . In particular, the images obtained with the inverse of 
Rivier's rule are not given by the average of the phase space representations 
using the inverses of the standard and antistandard rules. 

The arguments for demonstrating the feasibility of the inverse trans- 
formation F — > q are analogous, based on expanding F in power series 
J2nm b n m,q n p m and transforming each q n p m term independently to obtain the 
density operator as ^Z nm b nm Q nm - Also in this case the singularity is avoided 
due to the delta functions. The result is 



m n n—k 
1=0 k=0 j=0 

d k D, 



m 
I 



n 
k 



n — k 
J 



Hi 



l+ki 



■ l) l 2 k - n 



x 



d6 k 



0=0 



Qjpm-l-n-k-j 



(28) 



Moreover, using (|T|) and (|T6| ) it is easy to check that ( |T0D is verified. In 
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other words, Rivier's rule, as well as other rules based on functions / with 
derivatives at the origin have an inverse basis in a generalized sense. Eq. ([jjD 
is a symbolic expression that results from trying to fit the transformation (|5|), 
whose calculation is possible as shown above, into the scheme of eq. fll4|) . 
The fit requires a formal change in the order of integration. This generally 
illegitimate procedure is allowed when dealing with generalized functions. In 
this regard it is worth recalling that the integral form of the a delta func- 
tion [(27r) _1 J dxexp (ixy)} is in fact a symbol that is not to be interpreted 
"literally" . The origin of this symbol is again a formal change of integration 
when performing two successive Fourier transformations. In the actual com- 
putation the order of integration is reversed, see the illustrative discussion 
in |[25|| , or a rigorous and more general analysis in [^6| -especially sec. 7.9-. 
In the same vein, the actual order of integration when using (|15|) is to be 
reversed so that it is never performed over 9 and r first. This avoids the 
possible singularities of f^ 1 at the zeroes of /. 

In summary, the number of phase space formalisms equivalent to the con- 
ventional operator approach is broader than it had been generally believed. 
Zeroes of / do not preclude the possibility of a set of biunivocal transforma- 
tions between operators and phase space representations. 
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TABLE 1 



f 


F 


quantization rule 


1 


Wigner 


Weyl 


cos(9rh/2) 


Margenau-Hill 


Rivier(symmetrization) 


2sm(9Th/2)/9Th 


Shankara 


Born- Jordan 


e -i(0rh/2) 


Kirkwood 


standard 


e i(0rh/2) 


Kirkwood 


antistandard 




P-function(Sudarshan-Glauber) 


normal 




Q-function(Husimi) 


antinormal 



TABLE 1 CAPTION: A is a real parameter different from zero. 
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